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1 Introduction 



Nucleus-nucleus (A+A) collisions at relativistic energies have been intensely studied over 
the last two decades. The main goal of these efforts is to understand the properties of 
strongly interacting matter under extreme conditions of high energy and baryon densi- 
ties when the creation of the quark-gluon plasma (QGP) is expected [1]. Experimental 
results obtained in a broad range of collision energies indicate that an extended zone of 
strongly interacting dense matter indeed occurs in the collision process. Various collision 
characteristics and their collision energy dependence suggest [2, 3] that a transient state 
of deconfined matter may be created at collision energies as low as 40 AGeV. 

QGP formation is expected to occur at the early collision stage when the system is 
sufficiently hot and dense. In the course of further evolution, the system dilutes and 
cools down, hadronizes and finally decays into free hadrons. Therefore the final state 
hadrons carry only indirect information about the early stage of the collision. Thus firm 
conclusions about the creation of deconfined matter require studying a variety of collision 
characteristics. 

As fluctuations are sensitive to the dynamics of the system, in particular at the phase 
transition, the analysis of event-by-event fluctuations has been proposed as an important 
tool in the study of A-l-A collisions [4]. Large acceptance detectors, which allow the 
observation of a significant fraction of the final state particles, have made this suggestion 
truly attractive [5]. First pioneering studies were carried out by the NA49 collaboration 
on the fluctuations of the average transverse momentum [6] and the K/tt ratio [7] in 
central Pb+Pb collisions at the CERN SPS. 

There are a number of collision characteristics which can be inferred from experimen- 
tal data on event-by-event fluctuations. First of all, the fluctuation analysis can help to 
resolve the long standing problem whether, or rather to what extent, the strongly inter- 
acting matter, emerging from the early collision stage, achieves both thermal [8, 9, 10] 
and chemical [11, 12] equilibrium. In the NA49 investigation of event- wise fluctuations in 
the K/n ratio [7] in central Pb-|-Pb collisions at 158 AGeV, no signiflcant deviations from 
the hadro-chemical equilibrium ratio were found. If the equihbrium is indeed reached, the 
system's heat capacity [13, 14] and its compressibility [15] can, in principle, be deduced 
from the temperature and multiplicity fluctuations respectively. While large multiplicity 
fluctuations can be treated as a signal of particle production via clusters [16] or super- 
cooled droplets of deconfined matter [17], small final state fiuctuations of the conserved 
charges, electric or baryonic, can occur when fiuctuations generated in the plasma phase 
are frozen due to the rapid expansion of the system [18, 19]. On the other hand, significant 
transverse momentum and multiplicity fluctuations can result if the system hadronizes 
from a QGP near the predicted second order critical QCD end point [14, 20]. The latter 
has been shown by recent lattice QCD studies to occur at a substantial baryochemical 
potential [21], characteristic of the CERN SPS energy range. 

This paper extends the previous study [6] of the NA49 experiment of event-by-event 
transverse momentum fluctuations. The main objective is to observe how the fluctuation 
pattern changes with increasing number of nucleons participating in a collision, i.e. with 
the system size. For this purpose, not only Pb-|-Pb collisions at various centralities are 
studied, but also p-|-p, C-l-C and Si-|-Si interactions at a beam energy of 158 GeV per 
nucleon. In particular, it will be checked, whether the correlations present in the flnal 
state of p+p interactions survive in the collisions of heavier systems, as expected if the 
nucleus-|-nucleus (A-l-A) collision is a simple superposition of nucleon-|-nucleon (N-l-N) 
interactions. Moreover a test will be performed of the reasonable expectation that the 
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fluctuations become more similar to those of an equilibrated system when the number of 
participating nucleons increases. 

Measuring event-by-event fluctuations in A-l-A collisions one should consider the in- 
fluence of two trivial sources of fluctuations. The flrst one is caused by event-by-event 
fluctuations of the collision geometry and the second one by the flnite number of particles 
(statistical fluctuations). The dynamical fluctuations of interest have to be extracted from 
the noise caused by these trivial sources. 

In such a situation, a suitable choice of statistical tools for the study of event-by- 
event fluctuations is really important. In this work, mainly the fluctuation measure $ 
introduced in [8] will be employed. However, other fluctuation measures have also been 
proposed and studied, e.g. crpj.^dyn [22], AcTp^ [23] and Spj, [24], which can all be related. 
$ equals exactly zero when inter-particle correlations are absent. It also eliminates 'ge- 
ometrical' fluctuations due to the impact parameter variation. Thus, $ is 'deaf to the 
statistical noise and 'blind' to the collision centrality. The $ measure was already used 
in the previous experimental study [6] of pr fluctuations. It was also calculated within 
various models of nuclear collisions [25, 26, 27, 28, 29, 30]. However, in these model 
considerations, the effects of experimental acceptance, which significantly infiuence the 
results, were usually not taken into account. Consequently, most of these results cannot 
be compared to the data. 

For better understanding of the structure of the correlations contributing to two- 
particle transverse momentum correlations (as proposed in [23]) are studied as well. A 
preliminary analysis was presented in [31]. 

This paper is organized as follows: In Sec. 2 the statistical tools used in this anal- 
ysis are introduced and briefly discussed. The experimental set-up and procedures are 
presented in Sec. 3. Experimental effects such as detector acceptance and two-track reso- 
lution are discussed in Sec. 4. The results on the system dependence of the pt fluctuations 
are presented in Sec. 5. Discussion of the results and their comparison with theoretical 
models is given in Sec. 6. A summary closes the paper. 



There are numerous observables which can be used to measure pt fluctuations in high 
energy collisions. A natural one is the distribution of the average transverse momentum 
of the events deflned as 



where is the multiplicity of accepted particles in a given event and pTi is the transverse 
momentum of the i-th particle. The distribution of M{pt) is usually compared to the 
corresponding distribution obtained for 'mixed events' in which the particles are inde- 
pendent from each other and follow the experimental inclusive spectra (the multiplicity 
distribution for mixed events is the same as for the data). A difference between the two 
distributions signals the presence of dynamical fluctuations. Since the M{pt) distribu- 
tion crucially depends on the particle multiplicity, the method cannot be used to compare 
systems of signiflcantly different multiplicities. 

A more appropriate measure is the quantity $ [8] which, by its construction, is in- 
sensitive to the system size. In this paper the $ measure is used for the analysis of Pt 
fluctuations (^p^)- Following the authors of [8], one deflnes the single-particle variable 



2 Measures of fluctuations 




i=l 



(1) 



4 



Zpj, — Pt —Pr with the bar denoting averaging over the single-particle inclusive distribu- 
tion. One easily observes that = 0. Further, one introduces the event variable Zp^, 
which is a multi-particle analog of Zp^ , defined as 

N 

ZpT =J2iPTi-PT), (2) 
1=1 

where the summation runs over particles in a given event. Note, that (^pj,) = 0, where 
(...) represents averaging over events. Finally, the ^p^ measure is defined as 



(3) 



The second part of Eq. (3) is simply the dispersion of the inclusive Pt distribution 
(further denoted as cTpj,)- It can be easily shown that = 0, when no inter-particle 
correlations are present and the single-particle spectrum is independent of multiplicity. 
As already mentioned, is insensitive to centrality. This property may be expressed 
as follows: $pj, is independent of the distribution of the number of particle sources if the 
sources are identical and independent from each other [8, 10]. In particular, does not 
depend on the impact parameter if the A-|-A collision is a simple superposition of N-l-N 
interactions. 

In spite of the above mentioned advantages, there is an important disadvantage of 
using in the fluctuation analysis. While $pj, is sensitive to the presence of dynamical 
correlations in a system, it does not disentangle their nature. Fluctuations of very different 
character contribute to ^pj,. In order to achieve a better understanding of the fluctuation 
structure one needs to apply a more differential method [23]. 

The correlations can be studied by plotting the cumulative pr variables of particle 
pairs. Namely, for a given particle, instead of its px one introduces the variable x deflned 
as [32] 

fPT 

x{pt) = / p{pT')dpT', (4) 
Jo 

where p{pt) is the inclusive pr distribution, normalized to unity, which is obtained from 
all particles used in the analysis. By construction, the x variable varies between and 
1 with a flat probability distribution. The two-particle correlation plots, as presented in 
this paper, are obtained by plotting {xi,X2) points for all possible particle pairs within 
the same event. The number of pairs in each {xi, X2) bin is divided by the mean number 
of pairs in a bin (averaged over all {xi,X2) bins). This two-dimensional plot is uniform 
when no inter-particle correlations are present in the system. Correlations due to the Bose 
statistics produce a ridge along the diagonal of the {xi,X2) plot, which starts at (0, 0) and 
ends at (1, 1), whereas temperature fluctuations lead to a saddle shaped structure [23]. 

As will be seen in the flgures, the distribution of xi or X2 obtained from the two- 
dimensional {xi,X2) plots by projecting on the Xi or X2 axis is not flat. This is due to the 
method by which the plots are constructed. Namely, each pair of particles is represented 
by a point on the plot. Therefore, the events with higher multiplicities are represented 
by a larger number of pairs than those with smaller multiplicities. Since the shape of the 
Pt distribution depends on the event multiplicity, the projection of the two-dimensional 
plot on xi or X2 is no longer flat. However, it should be stressed that in the absence of 
any correlations the {xi,X2) plot is uniformly populated and the Xi and X2 projections 
are flat. 
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3 Experimental Set-up 



The NA49 experiment is a large acceptance liadron spectrometer at the CERN-SPS used 
to study the hadronic final states produced by collisions of various beam particles (p, Pb 
from the SPS and C, Si from the fragmentation of the primary Pb beam) with a variety of 
fixed targets. The main tracking devices are four large volume Time Projection Chambers 
(TPCs) (Fig. 1) which are capable of detecting 80% of approximately 1500 charged 
particles created in a central Pb+Pb collision at 158 AGeV. Two of them, the Vertex 
TPCs (VTPC-1 and VTPC-2), are located in the magnetic field of two super-conducting 
dipole magnets (1.5 and 1.1 T, respectively) and two others (MTPC-L and MTPC-R) 
are positioned downstream of the magnets symmetrically to the beam line. The results 
presented here arc analyzed with a global tracking scheme [33], which combines track 
segments that belong to the same physical particle but were detected in different TPCs. 
The NA49 TPCs allow precise measurements of particle momenta p with a resolution of 
a{p)/p'^ = (0.3 - 7) • 10-^ (GeV/c)-^ The set-up is supplemented by two Time of Flight 
(TOF) detector arrays and a set of calorimeters. 

The targets, C (561 mg/cm^). Si (1170 mg/cm^) and Pb (224 mg/cm^) for ion collisions 
and a liquid hydrogen cylinder (length 20 cm) for elementary interactions, are positioned 
about 80 cm upstream from VTPC-1. 

Pb beam particles are identified by means of their charge as seen by a Helium Gas- 
Cherenkov counter (S2') and p beam particles by a 2 mm scintillator (S2). Both of these 
are situated in front of the target. The study of C+C and Si+Si reactions is possible 
through the generation of a secondary fragmentation beam which is produced by a primary 
target (1 cm carbon) in the extracted Pb-beam. With the proper setting of the beam line 
magnets a large fraction oiall Z/A— 1/2 fragments at ^ 158 AGeV are transported to the 
NA49 experiment. On-line selection based on a pulse height measurement in a scintillator 
beam counter (S2) is used to select particles with Z = 6 (Carbon) and Z = 13, 14, 15 (Al, 
Si, P). In addition, a measurement of the energy loss in beam position detectors (BPD- 
1/2/3 in Fig. 1) allows for a further selection in the off-line analysis. These detectors 
consist of pairs of proportional chambers and are placed along the beam line. They also 
provide a precise measurement of the transverse positions of the incoming beam particles. 

For p, C and Si beams, interactions in the target are selected by an anti-coincidence of 
the incoming beam particle with a small scintillation counter (S4) placed at the beam axis 
between the two vertex magnets. For p+p interactions at 158 AGeV this counter selects 
a (trigger) cross section of 28.5 mb out of 31.6 mb of the total inelastic cross section. 
For Pb beams, an interaction trigger is provided by an anti-coincidence with a Helium 
Gas-Cherenkov counter (S3) directly behind the target. The S3 counter is used to select 
minimum bias collisions by requiring a reduction of the Cherenkov signal by a factor of 
about 6. Since the Cherenkov signal is proportional to Z^, this requirement ensures that 
the Pb projectile has interacted with a minimal constraint on the type of interaction. 
This setup limits the triggers on non-target interactions to rare beam-gas collisions, the 
fraction of which proved to be small after cuts, even in the case of peripheral Pb-|-Pb 
collisions. 

The centrahty of the nuclear collisions is selected by use of information from a down- 
stream calorimeter (VCAL), which measures the energy of the projectile spectator nucle- 
ons. The geometrical acceptance of the VCAL calorimeter is adjusted in order to cover 
the projectile spectator region by the setting of the collimator (COLL). 

Details of the NA49 detector set-up and performance of tracking software are described 
in [34]. 
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4 Data selection and analysis 



4.1 Data sets 

The data used for the analysis consists of samples of p+p, C+C, Si+Si and Pb+Pb 
collisions at 158 AGeV. For Pb+Pb interactions a minimum bias trigger was used allowing 
a study of the centrality dependence. The distribution of energy measured in the VCAL 
for the minimum bias Pb+Pb events was divided into six centrality bins (Table 1) which 
are numbered from 1 (the most central) to 6 (the most peripheral). For each bin of 
centrality the range of the impact parameter b, and the mean number of wounded nucleons 
(Nw) were determined by use of the Glauber model and the VENUS event generator [35]. 
The fraction of the total inelastic cross section of nucleus+nucleus collisions {(j/(Jtot) 
corresponding to each data set was calculated directly by use of the distribution of energy 
measured in the VCAL. In order to estimate the correlation between the energy deposited 
in the VCAL and the impact parameter minimum bias VENUS events were processed 
through the GEANT detector simulation code, and the energy deposited in the VCAL 
was simulated. The correlation between b and {Nw) was obtained from Glauber model 
calculations using the spectator-participant model of A+A interactions. The values of 
c/^tot; i^w) s-nd b presented in Table 1 are taken from [36] for minimum bias Pb+Pb 
collisions and from [37] for C+C and Si+Si interactions. 

The minimum bias Pb+Pb data consists of data taken during three different periods 
and at both magnetic field polarities. The most central Pb+Pb events correspond to 5% 
of the total geometric cross section. Since the minimum bias data provide only a small 
number of events in the most central Pb+Pb bin, additional central trigger runs were 
used. 

4.2 Event and particle selection 

The aim of the event selection criteria is to reduce a possible contamination with non- 
target collisions. The primary vertex was reconstructed by fitting the intersection point 
of the measured particle trajectories. Only events with a proper quality and position of 
the reconstructed vertex are accepted in this analysis. The vertex coordinate z along the 
beam has to satisfy [-2 — 2;o| < A.z, where the nominal vertex position zq and cut parameter 
Az values are: -579.5 and 5.5 cm, -579.5 and 1.5 cm, -579.5 and 0.8 cm, -578.9 and 0.4 
cm for p+p, C+C, Si+Si and Pb+Pb collisions, respectively. 

In order to reduce the contamination of particles from secondary interactions, weak 
decays and other sources of non- vertex tracks, several track cuts are applied. The accepted 
particles are required to have measured points in at least one of the Vertex TPCs. A cut 
on the so-called track impact parameter, the distance between the reconstructed main 
vertex and the track in the target plane, is applied {\bx\ < 2 cm and \by\ < 1 cm) to 
reduce the contribution of non-vertex particles. Moreover, particles are accepted only 
when the potential number of points (calculated on the basis of the geometry of the 
track) in the detector exceeded 30. The ratio of the number of points on a track to the 
potential number of points is required to be higher than 0.5 in order to avoid the counting 
of track segments instead of whole tracks. Only forward rapidity tracks (4.0 < < 5.5, 
rapidity calculated assuming the pion mass for all particles) with 0.005 < pr < 1-5 GeV/c 
are used in this analysis. 

The NA49 detector provides a large (but not complete) acceptance in the forward 
hemisphere. Two example plots of px versus azimuthal angle (see Fig. 1 for definition) 
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for positively charged particles (for the standard polarity of the magnetic field) arc shown 
in Fig 2. The solid lines represent a parametrization of the acceptance limits by the 
formula: 

where the values of parameters A, B and C depend on the rapidity interval as given in 
Table 2. These values apply to negatively charged particles as well, provided in Eq. 5 
is replaced by 0/|0|(18O — |0|). Only particles within the analytical curves are used in 
this analysis. This well defined acceptance is essential for later comparison of the results 
with models and other experiments. 

The NA49 detector is able to register particles produced in a significantly wider ra- 
pidity range covering almost the whole forward hemisphere. It would be interesting to 
determine the correlation measure in the mid-rapidity region. Preliminary results 
for the rapidity range 2.9 - 4.0 were reported [38]. However, the azimuthal acceptance 
in this region is more limited and systematic uncertainties affecting are not yet fully 
understood. 



4.3 Corrections and error estimates 

The statistical error on $pj, was estimated as follows. The whole sample of events was 
divided into 30 subsamplcs. The value of (^p^. was evaluated for each subsample and the 
dispersion {D) of the results was then calculated. The statistical error of was taken 
to be equal to D/-\/30. 

The event and track selection criteria reduce the possible systematic bias of the mea- 
sured values. In order to estimate the remaining systematic uncertainty, the values of 
cut parameters have been varied within a reasonable range and the systematic error has 
been estimated as a half of the difference between the highest and the lowest value. 
In addition, results obtained from the analysis of data taken at two different magnetic 
field polarities as well as from different running periods have been compared. 

Event cuts are used to reject possible contamination of non-target interactions, how- 
ever there is always a small fraction of remaining non-target events which can influence 
the $pj, values. The dependence of on the event selection cut is shown in Fig. 
3. The observed variation of with is small. The estimated systematic error is 
smaller than 0.85 MeV/c for peripheral Pb-|-Pb collisions, 0.55 MeV/c for Si-|-Si data and 
0.5 MeV/c for p-|-p events. 

The majority of tracks selected by the track selection criteria are main vertex tracks 
and the remaining fraction (ailO%) originates predominantly from weak decays and sec- 
ondary interactions with the material of the detector. In order to estimate the influence 
of this contamination on the measured value of ^pj., the impact parameter cut was var- 
ied (Fig. 4). A small increase of with increasing impact parameter cut is observed 
and may be due to the increased contribution of non- vertex tracks from weak decays and 
secondary interactions. By use of the VENUS [35] simulations these tracks were found to 
be correlated in pt thus increasing the measured value of ^p^ . The estimated systematic 
error due to the contamination of non- vertex tracks is smaller than 1.6 MeV/c for central 
Pb-|-Pb collisions, 0.75 McV/c for Si+Si data and 0.35 McV/c for p+p events. 

Losses of tracks due to the reconstruction inefficiency and track selection cuts influence 
the measured values. In order to estimate this effect, the dependence of on the 
percentage of randomly rejected particles was calculated. These dependences for the 
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most peripheral (6), the most central (1) Pb+Pb collisions and for p+p interactions are 
shown in Fig. 5. Within the considered kinematic region (forward rapidity) the tracking 
efficiency of our detector is higher than 95%. Fig. 5 implies that the bias due to tracking 
inefficiency is not higher than 0.5-1.0 MeV/c. 

As an estimate of the systematic error on a maximal error resulting from the above 
study has been taken. The systematic error is about 1.6 MeV/c for Pb+Pb collisions and 
1.2 MeV/c for p+p, C+C and Si+Si interactions. 

It has already been shown [6] that the limited two track resolution influences the 
measured values. In order to estimate this contribution several samples of mixed 
events (for different A+A collisions) were produced. Mixed events were constructed from 
original events, the multiplicities of mixed events being the same as in the case of real 
events but each particle in a mixed event taken at random from a different real event. 
The value calculated for the sample of mixed events was consistent with zero. In 
the second step the mixed events were processed by the NA49 simulation software. The 
resulting simulated raw data were reconstructed and the measure calculated. The 
obtained <I>py values are negative as expected for the anti-correlation introduced by the 
losses due to the limited two track resolution. The additive two track resolution correction 
is calculated as the difference (A^p^) between the values of ^p^ after and before this 
procedure. Fig. 6 presents this correction versus mean multiplicity. The lines correspond 
to an analytical parametrization of this dependence. The absolute values of the track 
resolution corrections are larger for heavier colliding systems where the density of tracks 
is relatively high. The absolute values of A^p^ are also larger for positively charged 
particles than for negatively charged ones, which is mainly due to higher track density for 
positive particles caused by signiflcantly larger number of protons than anti-protons. The 
A$py values are negative indicating that ^p^, measured with an ideal detector would be 
higher. For a given multiplicity, corrected for the limited two track resolution effect 
equals 'raw' minus the corresponding A^p^. 

5 Results 

The results shown in this section refer to accepted particles, i.e. particles that are accepted 
by the detector and pass all kinematic cuts and track selection criteria. The data cover 
a broad range in pt (0.005 < pr < 1-5 GeV/c). The rapidity of accepted particles is 
restricted to the interval 4.0 to 5.5 which corresponds to forward rapidities in the collision 
of equal mass nuclei (at 158 AGeV energy the center of mass rapidity equals 2.9 for flxed 
target geometry), where the azimuthal acceptance is large. 

The mean multiplicities of accepted particles, the dispersions = {N"^) — (N)^ of 
the multiplicity distributions, the mean inclusive transverse momenta, the dispersions ap^, 
of inclusive transverse momentum distributions and ^p^ values for all data sets used in 
this analysis are given in Table 3. The values shown in this table have been calculated 
for all accepted charged particles as well as for the negatively and the positively charged 
particles separately. All values of were corrected for the two track resolution effect. 

Fig. 7 shows the distributions of the mean (per event) transverse momentum M{pt) for 
p+p, Si+Si and central Pb+Pb collisions. Points correspond to data and the histograms 
to mixed events. The data are not corrected for any experimental effects. Events with zero 
accepted particle multiplicity are not taken into account. The small difference between 
M{pt) distributions for real and mixed events demonstrates that dynamical fluctuations 
are small. Moreover, no distinct class of events with unusual fluctuations is observed. 
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The width of the M{pt) distribution strongly decreases with the coUiding system size as 
expected from the increasing particle multiplicity. 

The fluctuation measure $pj, is more sensitive to small dynamical fluctuations. The 
measured values, corrected for two-track resolution, are plotted in Fig. 8 versus mean 
number of wounded nucleons for all accepted charged particles and also for positively 
and negatively charged particles separately. The values are small (when compared 
to (Tpy) for all investigated systems < 10 MeV/c), but a significant centrality 

dependence is observed. The values increase with the system size up to the maximum 
value which is observed for the most peripheral Pb+Pb coUisions, and then values 
decrease with increasing number of wounded nucleons. The value for Si+Si collisions 
seems to be lower than that for the most peripheral Pb+Pb ones although the number 
of wounded nucleons in both reactions is similar. This might suggest that {Nw) does 
not fully determine the fluctuations in A+A collisions. The measure calculated for 
positively charged particles is always lower than that for the negatively charged ones. The 
$pj, value for all charged particles is always higher than that for either the negatively or 
the positively charged particles. 

Two-particle correlation plots of the cumulant transverse momentum variable x are 
presented in Fig. 9 for p+p, C+C, Si+Si and three ccntralitics of Pb+Pb interactions 
(note the different color scales). It is seen that the plots arc not uniformly populated. In 
particular, significant long range correlations of about 40% (the color scale varies from 
0.75 to 1.6) are observed for the p+p data. This is rather unexpected when compared 
to the low value. These correlations are not seen when heavier colliding systems 
are studied. Instead, short range correlations become visible as an enhancement of the 
point density in the region close to the diagonal. They are most prominent for central 
Pb+Pb coUisions and are consistent with the effect of Bose-Einstein statistics. For the 
most peripheral collisions the pattern seen in the two-particle correlation plot is different 
from that in the remaining systems. One observes an enhancement in the region close to 
xi = 1 and X2 = 1. 

6 Discussion 

In this section the results are compared with predictions of models and with the results 
of other experiments. 

Fig. 8 shows that $pj, is a non-monotonic function of centrality with the maximum 
at approximately N]^ = 40. Such a behavior strongly resembles the dependence of the 
magnitude of collective fiow - directed (vi) and elliptic (^2) - on Nw [39]. So, there 
is a natural suggestion that the px fiuctuations measured by may be caused by 
the collective fiow. This suggestion was checked by performing a simple Monte Carlo 
analysis generating events with independent particles, following the measured inclusive 
Pt distribution. The uniform azimuthal angle distribution of the events was modified by 
the collective elliptic fiow but the particles remained independent from each other with 
respect to their transverse momenta. For such events the dynamical px fiuctuations vanish 
if the azimuthal angle acceptance is complete. However, the azimuthal angle acceptance 
of the NA49 detector is not fiat, and consequently the azimuthal anisotropy generates 
a finite value of ^pj,, even though the px of particles are independent of each other. 
It was found that an abnormally large value of ^2 = 0.5 leads to <l>py as large as 17 
MeV/c, but a realistic magnitude of V2, which depends on px (according to [39] taken as 
0.05 •pT(GeV/c)), results in a value of which is consistent with zero ($pj, = 0.0 + 0.6 
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MeV/c). Thus, one concludes that the effect of the azimuthal anisotropy caused by the 
collective flow combined with the incomplete azimuthal acceptance is not responsible for 
the observed dynamical pt fluctuations. 

In Fig. 10 the dependence of on the mean number of wounded nucleons is directly 
compared to predictions of the HIJING [40] model (default parameters were used) for all 
charged particles, and for negatively charged and positively charged particles separately. 
The same kinematic cuts are applied as for the data. The black lines represent the results 
of the HIJING simulations where the effect of the limited NA49 acceptance (pr versus 
azimuthal angle) is taken into account. The gray hues refer to the HIJING predictions 
for full azimuthal acceptance to demonstrate the effect of the hmited acceptance of the 
detector. 

In contrast to the data, $pj, computed within the HIJING model does not change when 
going from elementary to central Pb+Pb collisions because the HIJING model represents 
an essentially independent superposition of N+N interactions. The effects of short range 
correlations (Bose- Einstein and Coulomb) have not been incorporated in the HIJING 
model. However it was estimated in the previous analysis [6] that the combined effect of 
short range correlations produces values on the level of 5 MeV/c for central Pb+Pb 
collisions. This effect strongly depends on multiplicity and becomes negligible for p+p 
interactions. 

For the HIJING model values for positively charged particles are, as in the case 
of real data, lower than for negatively charged and for all charged particles. The fact that 
values for positively charged particles are always lower than for negatively charged 
ones has been found (using HIJING) to be related to the limited acceptance and treating 
protons as pions (assuming the pion mass for all produced particles). 

When was first introduced [8] it was believed that its value would be non-zero for 
elementary interactions (mainly due to the M{pt) versus dependence) and would vanish 
for heavier colliding systems as a result of equilibration. The present measurements do 
not confirm this expectation. Although is close to zero for central Pb+Pb collisions, 
the maximum value is observed not for p+p data, but for colliding systems with Nw — 
40. 

Although the value of $pj, is small for p+p collisions, a significant structure appears 
in the two-particle correlation plot (Fig. 9). The first candidate for its origin is the 
dependence of M{pt) on A^, observed for elementary interactions [8]. Fig. 11 compares 
{M{pt)) versus A^ for the HIJING model and for real p+p data ((...) represents averaging 
over events with a given A^). The HIJING model includes the NA49 acceptance and 
all kinematic restrictions. It shows good agreement with the measurements. Fig. 12 
presents the p+p two-particle correlation plots for data (a), for the HIJING model (b) 
and for a simple random generator, which reproduces the dependence of M{pt) on A^ 
observed in the data (c). The accepted particle multiplicity distribution for the random 
generator (c) is the same as in the case of data (all kinematic cuts and NA49 geometric 
acceptance are included). Both models (b) and (c) qualitatively reproduce the structure 
of the two-particle correlation plot observed in the data, however, the HIJING model (b) 
shows additionally a small enhancement of the point density in the region of high x. 

The value calculated for model (c) for all charged particles equals 1.2 ± 0.2 ± 1.8 
MeV/c and is consistent with for the p+p data (2.2 ± 0.3 ± 1.2 MeV/c). Nevertheless 
the small difference might indicate that, in agreement with the analysis presented in [29], 
there could be an additional source of correlations present in the data. However, due 
to the relatively high statistical and systematic errors the effect cannot be estimated 
quantitatively. A similar analysis was carried out for negatively charged particles only. 
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where the correlations caused by resonances and by charge conservation are expected to 
be smaller than for all charged particles. There, the experimental value of is 0.8 ± 0.1 
± 1.2 MeV/c while the model (c) gives 0.6 ± 0.2 ± 0.3 MeV/c. Two-particle correlation 
plots for negatively charged particles only are very similar for data and for model (c). 
Thus, one concludes that the results on negative particles from p+p interactions are 
consistent with the conjecture [41] that the particles are emitted independently but that 
their y and px distributions depend on event multiplicity. 

The last panel (d) of Fig. 12 presents the result of a simple temperature fluctuation 
model (the concept is described in [30]), which assumes that the only source of fluctuations 
is event-by-event fluctuation of the inverse slope parameter (T) of the transverse mass 
spectra. The model assumes a Gaussian shaped rapidity distribution and an exponential 
shape of the transverse mass distribution with the mean inverse slope parameter (T) = 
152 MeV adjusted to agree with the experimental p-|-p results. All kinematic cuts applied 
for the real data are also used in this model and the effect of the finite detector acceptance 
is taken into account. The mean multiplicity of all accepted particles is the same as in the 
data. The fluctuations of the inverse slope parameter lead to a saddle shaped structure 
in the two-particle correlation plots. Panel d) shows the result for fluctuations of T on 
the level of about 10% (the dispersion ctt = 16 MeV). Because of the difference between 
the panels a) and d) one concludes that ffuctuations of the inverse slope parameter are 
not the (main) source of correlations in p-|-p data. 

The HIJING model has also been used to obtain a two-particle correlation plot for 
G+G collisions (Fig. 13 b), which appears to be similar to that observed for real events 
(Fig. 9). Fig. 13 a (the same as Fig. 12 b) presents p+p events simulated by the HlJlNG 
model. The structure observed for p+p collisions vanishes for heavier systems due to the 
dilution effect from the higher number of uncorrelated particles (resulting from different 
N+N interactions), whereas the $pj, measure is not affected by this dilution effect. 

In order to see how dynamical fluctuations influence two-particle correlation plots for 
central Pb+Pb data the above model with fluctuations of the inverse slope parameter 
was used again. The mean inverse slope parameter (T) was set to 190 MeV and the 
mean multiplicity of all accepted particles to 200 in order to compare the results with 
central Pb+Pb collisions. The inverse slope parameter varied from event to event with 
a Gaussian shaped distribution of width a-r- Fig. 14 presents {xi,X2) plots for different 
levels of the inverse slope parameter fluctuations. The fluctuations lead to a saddle shaped 
structure which is not visible in central Pb+Pb collisions. One can thus exclude signiflcant 
T fluctuations in central Pb+Pb collisions at top SPS energy. In Fig. 15 the predicted 
dependence of on T fluctuations [30] is plotted and compared to measured for 
the 5% most central Pb+Pb interactions. The experimental value contains both 
short and long range correlations. The sohd fine corresponds to the T-fluctuation model 
presented in [30] which does not include short range correlations. The dashed line is 
the combination of this model with a contribution of short range correlations estimated 
experimentally [6] as 5 MeV/c, in agreement with theoretical arguments [9, 10]. One sees 
that the observed value of is already below the contribution of the Bose-Einstein 
correlations and that the inclusion of slope fluctuations makes the difference even larger. 
Thus, one can conclude, in agreement with the previous results [6], that the data leave 
no space for signiflcant T fluctuations provided they are not canceled by other negative 
correlations. 
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An increase of transverse momentum fluctuations was predicted [14] to occur in A+A 
collisions which freeze out near a second order critical end point of the QCD phase di- 
agram. Based on calculations from [14] and the numbers given in Table 3 it can be 
estimated that such critical fluctuations alone should result in ~ 20 MeV/c. This 
number is significantly larger than the maximum value of found in this analysis - for 
peripheral Pb+Pb interactions ($pj, = 7.2 ± 0.7 ± 1.6 MeV/c). Note, however, that in 
this theoretical estimate the effect of the limited experimental acceptance was not taken 
into account. 

Transverse momentum fluctuations in A+A collisions and elementary interactions were 
measured by several experiments in the SPS energy range. A value of = 10.9 ± 1.5 
MeV/c for charged hadrons was reported by the NA22 experiment in n^p and K~^p 
interactions [42] at 250 GeV in a rapidity acceptance similar to the one used in this 
paper. The value obtained by the NA22 experiment is higher than the result presented 
in this paper for p+p interactions at 158 AGeV = 2.2 ± 0.3 ± 1.2 MeV/c). This 
can be caused by the differences in azimuthal acceptance, energy and types of interacting 
hadrons. 

The CERES experiment at the SPS measured 5 MeV/c for charged hadrons in 

central Pb+Au collisions at 158 AGeV [43]. This measurement was performed close to 
mid-rapidity (pseudo-rapidity 2.2 < rj < 2.7) in the pr range 0.1 < pr < 1-5 GeV/c and 
with full azimuthal angle acceptance. A somewhat smaller value, 1.5 MeV/c, is 

found in this paper for the similar Pb-|-Pb reaction at forward rapidities. 

Signiflcant non-statistical fluctuations {^pj. = 52.6 ± 0.3 MeV/c) were measured for 
charged hadrons at mid-rapidity in central Au-|-Au collisions at ^/s^n = 130 GeV by the 
STAR experiment [44]. At the RHIC energies Pt fluctuations arc expected to be larger 
than at the SPS due to a signiflcant contribution of correlated particles originating from 
(mini-)jet fragmentation. 

A non-monotonic centrality dependence of pt fluctuations, with a maximum for semi- 
central Au-|-Au collisions, was reported by the PHENIX experiment at ^/s^n — 200 GeV 
[45]. This result is in qualitative agreement with the result presented in this paper for 
Pb+Pb collisions at 158 AGeV. A centrality dependence of pt fluctuations was observed 
also by the STAR experiment at ^/sNN = 130 GeV [44]. 

It should be stressed that currently available event-by- event measurements were per- 
formed in a limited acceptance, the multiplicity of observed particles being below 20% of 
the total multiplicity. Consequently, the sensitivity of these measurements, in particular 
to long range correlations is reduced. It is desirable to perform future event-by-event 
measurements in an extended acceptance. 

7 Summary 

Transverse momentum cvcnt-by-event fluctuations were studied for p+p, C+C, Si+Si and 
Pb+Pb collisions at 158 AGeV. The analysis was limited to the forward rapidity region. 
Three different characteristics were measured: the ffuctuations of average transverse mo- 
mentum {M{pt)) of the event, the fluctuation measure, and transverse momentum 
two-particle correlations. All measured $pj, values are below 10 MeV/c and are much 
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smaller than the dispersion of the inclusive pt distribution. However, the correlations 
observed in p+p collisions are not simply more and more diluted when going to heavier 
colliding systems as could be expected if the created matter approaches a higher level of 
equilibrium with increasing system size. Instead, a significant system size dependence of 
the measure is seen with a maximum for peripheral Pb+Pb collisions with ~ 
40. The two-particle correlation plot for p+p data shows a prominent structure which 
was found to be connected with the dependence of M{pt) on A^. This structure disap- 
pears when going to heavier colliding systems. Instead short range correlations become 
visible as an enhancement of the point density in the region close to the diagonal. This 
effect is strongest for the most central Pb-|-Pb interactions. No structure characteristic of 
event-by-event temperature fluctuations is observed. 

The HIJING model qualitatively reproduces the structure of two-particle correlation 
plots for p+p and C+C data. However, in contrast to the data, it shows no centrality 
dependence of $pj,. 

In future, a study of the energy dependence of transverse momentum fluctuations in 
the CERN SPS energy range is planned using the NA49 Pb+Pb collision data taken at 
different beam energies. The aim is to search for possible anomalies connected with the 
onset of the deconfinement phase transition, which is indicated by features of pion and 
strangeness production at low SPS energies [3]. 
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No. of events 


c/c"tot in each bin 


{Nw) 


b range [fm] 


D+D 


570 000 


0.9 


2 




c+c 


33 000 


0.153 


14 


- 2.0 


Si+Si 


63 000 


0.122 


37 


- 2.6 


Pb+Pb(6) 


117 000 


0.57 


42 


10.2 - 


Pb+Pb(5) 


59 000 


0.10 


88 


9.1 - 10.2 


Pb+Pb(4) 


68 000 


0.10 


134 


7.4 - 9.1 


Pb+Pb(3) 


68 000 


0.11 


204 


5.3 - 7.4 


Pb+Pb(2) 


45 000 


0.075 


281 


3.4 - 5.3 


Pb+Pb(l) 


180 000 


0.05 


352 


- 3. 1 



Table 1: Data sets used in analysis. Listed for p+p, C+C, Si+Si and six centralities of 
Pb+Pb collisions at 158 AGcV arc: number of events, cr/atot - the fraction of the total 
inelastic cross section in that bin, {Nw) - the mean number of wounded nucleons, h - the 
impact parameter range. 



Vtt 




1- c J 




3.9 - 4.1 





0.3 


6500 


4.1 - 4.3 





0.3 


5500 


4.3 - 4.5 





0.25 


4500 


4.5 - 4.7 





0.25 


3500 


4.7 - 4.9 





0.2 


2500 


4.9 - 5.1 


0.5 


0.2 


2500 


5.1 - 5.3 


1.0 


0.1 


2500 


5.3 - 5.5 


1.5 


0.1 


2500 



Table 2: The parametrization of the NA49 y —pr acceptance at 158 AGeV for positively 
charged particles (standard configuration of magnetic field) . For negatively charged parti- 
cles one has to redefine the azimuthal angle (see text) and then use the same parametriza- 
tion. 
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(N) 




[MeV/c] 


cTpy [MeV/c] 


[MeV/c] 


p+p (all) 


1.4 


1.3 


304 


196 


2.2 ± 0.3 


p+p (-) 


0.6 


0.7 


283 


179 


0.8 ± 0.1 


P+P (+) 


0.8 


0.9 


317 


206 


-1.4 ± 0.3 


C+C (mU) 


10 


1.3 


300 


210 


5.1 ± 0.7 


C+C (-) 


4.5 


2.4 


279 


190 


1.8 ± 0.8 


C+C (+) 


5.5 


2.7 


317 


224 


0.7 ± 0.7 


Si+Si (all) 


27 


7 


301 


217 


4.9 ± 0.8 


Si+Si (-) 


12 


4 


277 


195 


2.6 ± 0.5 


Si+Si (+) 


15 


4 


320 


231 


-0.2 ± 0.7 


Pb+Pb(6) (all) 


39 


18 


299 


220 


7.2 ± 0.7 


Pb+Pb(6) (-) 


18 


9 


270 


195 


4.5 ± 0.5 


Pb+Pb(6) (+) 


21 


10 


325 


237 


1.9 ± 0.7 


Pb+Pb(5) (all) 


73 


17 


305 


226 


6.6 ± 0.7 


Pb+Pb(5) (-) 


34 


9 


273 


199 


4.5 ± 0.7 


Pb+Pb(5) (+) 


39 


9 


333 


245 


0.6 ± 0.8 


Pb+Pb(4) (all) 


104 


19 


309 


230 


5.6 ± 0.8 


Pb+Pb(4) (-) 


49 


10 


276 


202 


3.8 ± 0.5 


Pb+Pb(4) (+) 


55 


11 


337 


249 


-0.6 ± 0.9 


Pb+Pb(3) (all) 


148 


21 


312 


233 


4.6 ± 0.8 


Pb+Pb(3) (-) 


69 


11 


279 


204 


2.9 ± 0.8 


Pb+Pb(3) (+) 


79 


12 


342 


252 


-1.3 ± 0.8 


Pb+Pb(2) (all) 


193 


21 


315 


234 


2.2 ± 1.0 


Pb+Pb(2) (-) 


90 


11 


281 


205 


2.4 ± 0.8 


Pb+Pb(2) (+) 


103 


13 


344 


254 


-3.7 ± 1.1 


Pb+Pb(l) (all) 


230 


19 


317 


236 


1.4 ± 0.8 


Pb+Pb(l) (-) 


108 


11 


281 


203 


0.9 ± 0.6 


Pb+Pb(l) (+) 


122 


12 


349 


257 


-2.9 ± 0.8 



Table 3: Measured inclusive and event-by-event parameters for accepted particles. (A^), 
(Tjv, Pt and (Tp^, values are not corrected for acceptance. values arc corrected for 

limited two track resolution. The systematic error of is smaller than 1.6 MeV/c. 
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Figure 2: NA49 (p — pr acceptance of positively charged particles (standard configuration 
of magnetic field) for two selected rapidity bins at 158 AGeV. The solid lines represent 
the analytical parametrization. 
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Figure 3: The dependence of on the allowed distance Az of the reconstructed event 
vertex from its nominal position. Note: the values and their errors are correlated. 
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Figure 4: The dependence of on the upper cut in the impact parameter For each 
point the cut on \by\ was equal half the cut on \bx\. Note: the values and their errors are 
correlated. The dashed line indicates the cut used in the analysis. 
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Figure 5: The dependence of on the fraction of randomly rejected particles. 
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Figure 6: Additive correction A$py for limited two track resolution effect versus multiplic- 
ity (N) of accepted particles. Different points correspond to positively charged, negatively 
charged and all charged particles. The presented corrections have been determined us- 
ing p-|-p, C-l-C, Si-|-Si and three centralities of Pb-|-Pb coUisions at 158 AGeV. The hues 
represent the analytical parametrization: A^p^{{N)) — —a^ (N) + b with a and b being 
parameters of a fit to the data points. 
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Figure 7: Distributions of mean transverse momentum for real (data points) and mixed 
events (histograms). Data points are not corrected for acceptance and limited two track 
resolution. Events with accepted particle multiplicity equal to zero are not used. 
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Figure 8: versus mean number of wounded nucleons (Nw)- Data points were corrected 
for limited two track resolution. Errors are statistical only. Systematic error is smaller 
than 1.6 MeV/c. 
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Figure 9: Two-particle correlation plots using the cumulant pt variable x. After each 
charged particle pair {xi,X2) was entered into the plot, the bin contents were normalized 
by dividing with the average number of entries per bin. The data are plotted with different 
color scales. 
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Figure 10: versus mean number of wounded nucleons calculated using the HIJING 
model with geometrical acceptance cuts included (black lines) and without geometrical 

acceptance restrictions (gray lines). Results are compared to data (points) corrected 
for limited two track resolution (the markers are the same as in Fig. 8). The panels 
represent: all charged, negatively charged and positively charged particles. Data points 
contain both short and long range correlations. The effects of short range correlations are 
not incorporated in the HIJING model. 
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Figure 11: Mean transverse momentum versus multiplicity of all accepted particles. The 
closed symbols represent p+p data at 158 AGeV (data are not corrected for limited two 
track resolution effect) and the open symbols corresponds to p+p events simulated using 
the HIJING model (the effects of the limited NA49 acceptance are included). Events with 
accepted particle multiplicity equal zero are not used. 
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Figure 12: Two-particle correlation plots using the cumulant pt variable x. Results are 
shown for all charged particles for p+p data (a) compared to: simulated p+p HIJING 
events with limited NA49 acceptance (b), simple random-generator model, which repro- 
duces M{pt) versus correlation for p-t-p data (c), model of fluctuations of the inverse 
slope parameter for p+p data on the level of about 10 % (d). 
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Figure 13: Two-particle correlation plots using the cumulant px variable x. Results are 
shown for simulated p+p (a) and C+C (b) collisions from the HIJING model. Limited 
NA49 acceptance is taken into account. 
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Figure 14: Two-particle correlation plots using the cumulant px variable x for the most 
central Pb+Pb collisions. The experimental result is shown in (a) and compared to the 
inverse slope parameter fluctuation model (b-d). In the model the mean value of the 
inverse slope parameter was set to 190 MeV. The dispersions of the Gaussian shaped 
inverse slope parameter distributions were set to: ctt = 5 MeV (b), ar = 10 MeV (c), ax 
= 20 MeV (d) 
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Figure 15: Predicted dependence of ^pj, on T fluctuations [30] compared to the measured 
$pj, for the 5% most central Pb+Pb colhsions. The uncertainty of the measured 
value is represented by the bands for statistical (gray) and systematic (dark gray) errors. 
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